FINAL  REPORT  ON 
AFOSR  CONTRACT  F49620-83-C-0064 


Steven  A.  Orszag,  Principal  Investigator 
Department  of  Mathematics 
MIT 

Cambridge,  MA  02139 
Volume  1 


DTIC 


ELECTE 
MAR  0  11988  = 


Ok 


H 


PBFTRIBUTION  8TATEMPIT  A 

Approved  for  pid^Uc  leleoee; 
Dietrtbution  Unlimited 


3 


1 


UNCLASSIFIED 

ieCU«iTy  classification  of  th«s  page 


lA  «6PO«r  SECUfliTv  classification 

Unclassified 


2*.  SECLifliTV  CLASSIFICATION  AUTHORITY 


2b.  oeclassificatiqn/oowngraoing  scheoule 


4.  PERFORMING  organization  REPORT  NUMaERlSI 


REPORT  DOCUMENTATION  PAGE 


10.  RESTRICTIVE  MARKINGS 


3.  oistri0ution/availa8il'Tv  op  report 

-Qr  r'-4.l<ASf  ^ 

clidr]  »■  I  MH /I'fcM 


S.  monitoring  organization  report  NUMeER(S) 

AFOSR-TK-  8  7  -  1  ^  90 


6m.  NAME  OF  performing  ORGANIZATION  Bo.  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 

ilf  appiiemtttmp  , 


6e.  AOORE5S  (Cily.  Stmtm  and  ZIP  Codai 

Deoartment  of  t^lathenatics 
tllT 

Cambridae,  MA  02139 


•a.  NAME  0^  FUNOINC/SPONSaniNa 
OnGANIZATION 

afosr/na 


7b.  AOORESS  iCity.  Statm  and  ZIP  CoOmt 

Su.jUm<x  J  AF3  Od, 

^0  T.^2.  -  t  “ 


AFOSl.  /m A 


IBb.  OFFICE  SYMeOL  9.  PROCUREMENT  INSTRUMENT  rOENYlFlCATION  NUMBER 

(If  mppUeabUi  ,  , 

NA  F4S620-83-C-0064 


10.  SOURCE  OP  PUNOING  NOS. 


PROGRAM 

PROJECT 

TASK 

element  no. 

NO. 

NO. 

6//D5-  f 

A  2. 

13.  PEASONAu  auTmO«(S1 

-  Steven  A.  Orszaa 


Ua  Tvpe  op  REpo«T  13b.  time  covebeo 

Final  Report  pbom  2/1/83  to 


16.  supplementary  notation 


14.  date  of  report  fYr.,  Mo.,  Doyf  I  19.  PAGE  COUNT 


4  May,  1987 


COSATI  COOES 


16.  SUBJECT  TERMS  rConnnua  on  rwmrtm  if  nawiaoo'  •'Vd  idanrl^y  by  btoek  numbmrt 

”T”tAr  la K  I-*  j  (  eal  STim  k  la"!  |  oia.  '' 


19.  ABSTRACT  (Conttnuo  on  rmvonm  \f  nmttmmory  and  idanct/y  by  block  numbcrt 

We  sumnarize  work  done  under  AFOSR  Contract  F49620-83-C-0064.  TTie  major  results  include: 
development  of  renormalization  group  techniques  for  large-eddy  simulations  of  turbulent 
flews,  the  first  direct  numerical  simulation  of  turbulent  spots  in  channel  and  boundary- 
layer  flows,  the  further  development  of  spectral  methrxls  for  turbulence  simulations, 
the  identification  of  secondary  instability  modes  in  free  shear  layers,  the  development 
of  an  efficient  multi-grid  marching  method  for  solution  of  the  parabolized  Navier-Stokes 
equations,  a  mathematical  analysis  of  boundary  conditions  for  the  parabolized 
compressible  Navier-Stokes  equations,  the  further  development  of  a  method  to  improve 
numerical  solution  of  singular  perturbation  problems  by  use  of  asymptotic  approximations. 


I 

I 


m 


30.  QISTRI  auTION/AvAICAaiLlTV  op  AaSTRACT 
UNCI.AS3IPIEO/UNV.IMITEO  H  SAME  AS  RPT.  X  OTIC  USERS  H 


23a  name  op  RESPONSiauE  INOIVIOUAC 

Dr  vJa.MA*s  M  )^c.H'>cl\a.-c( 


00  FORM  1473,  83  APR  eo.t.cn 


21.  abstract  SECURITY  CLASSIFICATION 
A  C  /a  S  i  I  €cl 


22b.  TELEPHONE  NUMBER 
(inctudb  Arwo  Codbf 


22e.  OFFICE  SYMBOL 


r/nciud#  Aryo  Qodbt  t 

(ztf  I'N  74  7'4'‘7S  6  APdsA/u4 


EOiTiCN  OF  I  ,AN  73  IS  OBSOLETE. 


WWW 


FINAL  REPORT  ON  AFOSR  CONTRACT  F49620-83-C-0064 


Steven  A.  Orszag,  Principal  Investigator 
Department  of  Mathematics 

MIT 

Cambridge,  MA  02139 

‘  r 

'V 

— In  the -atiaehc4-papers.-^e  summarize  work  done  on  this  research 

project.  The  major  results  include: 

1.  Development  of  renormalization  group  techniques  for  large-eddy 
simulations  of  turbulent  flows. 

1,  The  first  direct  numerical  simulation  of  turbulent  spots  in  channel 
and  boundary  layer  flows. 

3.  The  further  development  of  spectral  methods  for  turbulence 
simulations. 

A.  The  identification  of  secondary  instability  modes  in  free  shear 
layers. 

5.  The  development  of  an  efficient  multi-grid  marching  method  for 
solution  of  the  parabolized  Navier-Stokes  equations. 

6.  A  mathematical  analysis  of  boundary  conditions  for  the  parabol¬ 
ized  compressible  Navier-Stokes  equations.  ^ 

7.  The  further  development  of  a  method  to  improve  numerical  solu¬ 
tion  of  singular  perturbation  problems  by  use  of  asymptotic  ap¬ 
proximations. 

Further  details  are  given  in  the  attached  papers. 
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SUMMARY 

The  Booster  Method  for  improvement  of  the  numerical  solution  of  par¬ 
tial  differential  equations  by  the  addition  of  asymptotic  corrections  to 
the  right  hand  side  is  presented.  It  is  applied  here  to  the  diffusion- 
convection  equation  for  the  case  of  'small'  diffusion.  The  correction 
terms  were  used  in  finite  difference  and  finite  element  schemes.  The 
finite  element  results  were  used  as  reference  for  checking  the  performance 
of  the  finite  difference  schemes.  Excellent  results  were  obtained  without 
the  use  of  upstreaming  or  artificial  diffusion.  Theoretical  expectations 
were  confirmed. 

1.  INTRODUCTION 

Singularly  perturbed  initial  and  boundary  value  problems  for  partial 
differential  equations  appear  in  various  fields  of  application  such  as 
fluid  dynamics,  heat  transfer,  transport  of  atmospheric  pollution,  etc.  In 
particular,  such  equations  appear  in  diffusion-convection  processes.  Often 
the  (normalized)  diffusion  coefficient  e  becomes  small,  and  thin  boundary 
or  interior  layers  appear  within  the  region  of  interest.  Consequently, 
these  problems  become  increasingly  difficult  to  solve  numerically  by  dis¬ 
cretization  methods. 

We  would  like  to  avoid  the  use  of  a  prohibitively  large  number  of  grid 
points,  as  required  for  resolution  by  straightforward  numerical  methods 
when  c  decreases.  To  this  end,  several  approaches  are  possible,  such  as 
the  use  of  nonuniform  meshes,  adaptive  techniques,  positive  type  schemes, 
etc.  The  question  of  applicability  of  such  schemes  to  multi-dimensional 
problems  is  presently  open. 

A  different  approach  is  motivated  by  classical  singular  perturbation 
methods  where  'inner'  and  'outer'  solutions  are  combined  to  give  approxi¬ 
mate  solutions.  These  solutions  become  more  accurate  as  the  equations 
become  stiffer,  however,  the  error  is  fixed  for  a  given  G  and  cannot  be 
improved  or  estimated  reliably  in  most  cases  of  interest. 

The  Booster  Method  attempts  to  combine  the  asymptotic  approach,  with 
known  discretization  methods,  in  order  to  obtain  a  numerical  method  which 
improves  when  e  becomes  smaller.  At  the  same  time,  it  keeps  the  proper¬ 
ty  that  the  error  can  be  made  arbitrarily  small  for  any  fixed  e  by  refin¬ 
ing  the  computational  mesh  (Israeli  and  Ungarish  [1],[2]).  For  the  one¬ 
dimensional  case,  we  were  able  to  prove  that  an  improvement  by  a  factor  of 
can  be  obtained  where  e  is  the  'small'  parameter  and  n  is  the 
order  of  the  asymptotic  approximation  used  (Israeli  and  Ungarish  [2])  .  We 
expect  similar  behaviour  in  the  multidimensional  case  [1]. 

In  the  present  paper,  we  investigate  a  multi-dimensional  application 
to  diffusion-convection  problems. 
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2 .  FORMULATION 

We  consider  the  transport  of  a  quantity  q  in  a  rectangular  region. 
The  normalized  partial  differential  equation  is 

L(q)  i  -EV^q  +  V*7q  =  0  .  (1) 

In  the  present  application  the  velocity  field  V  is  assumed  to  be 
known  and  q  is  specified  on  the  boundaries.  This  problem  was  often  used 
as  a  test  case  for  various  finite  difference  and  finite  element  methods  of 
solution  and  it  is  well  known  that  most  methods  fail  as  the  cell  Reynolds 
number  (lv|h)/e  becomes  larger  than  0(1)  (here  h  is  a  representative 
mesh  size) .  For  example  centered  schemes  develop  unphysical  oscillations 
in  space,  while  uncentered  schemes  have  unaccepta.ble  artificial  diffusion 
and  are  of  lower  order  over  the  same  computational  stencil. 

The  Booster  Method  uses  an  asymptotic  approximation  q(x,y)  to  the 
solution  q(x,y)  in  order  to  reduce  the  truncation  error  in  the  numerical 
scheme . 

The  'usual'  numerical  solution  Q(x,y)  (defined  only  at  grid  points) 
is  obtained  from 

Lj^(Q)  =  f  .  !2) 

where  Lfj  is  the  discrete  approximation  to  the  differential  operator  L. 
The  improved  numerical  solution  Q  is  obtained  from 

"  f  L  (q)  ~  I<(q)  .  (3) 

N  N 

Here  I,(q)  is  the  differential  operator  applied  to  the  approximate  solu¬ 
tion.  Thus  the  Booster  Method  applies  an  asymptotic  correction  to  the 
right  hand  side  of  the  equation  and  therefore  requires  a  negligible  amount 
of  extra  work.  It  can  be  used  with  any  numerical  scheme  without  modifica¬ 
tion  in  the  method  ojf  solution. 

The  same  basic  approach  of  using  asymptotic  corrections  to  the  right 
hand  side  can  be  used  to  improve  the  Standard  Finite  Element  (SFE)  method. 
The  resulting  Asymptotic  Finite  Element  (AFE)  method  is  described  briefly 
in  the  following;  for  details  see  Bar-Yoseph  and  Israeli  [4],  [5]. 

Suppose  that  the  unit  square  is  divided  into  elements  and  that  the 
variation  of  q  within  the  given  region  is  approximated  by 


Q(x)  »  N.  (x)Q.  ,  (4) 

i  1 

where  is  the  value  of  the  approximate  solution  at  the  i-th  nodal  point 

and  is  the  corresponding  global  trial  function  (we  use  the  summation 

convention,  with  summation  over  the  nodes  within  the  given  region) .  The 
Bubnov-Galerkin  finite  element  scheme  of  eq.  (1)  is  given  by 


e(7  n.,7n.Q.)^+  (N.  ,V*Vn,Q.)^ 

j  i  1  h  j  11  h 


(N.,f)^ 
3  h 


1,2, ... ,m 


where  m  is  the  number  of  inner  nodal  points  and  (*,*)  denotes  the  usual 
inner  product  in  1.2^^)  •  The  subscript  h  in  ‘ienotes  an  approxi¬ 
mation  to  (*,•)  obtained  by  a  quadrature  rule. 

Our  corresponding  asymptotic  finite  element  (AFE)  scheme  for  eq.  (1) 
is  the  following 


G(V  N.,VN.Q.)^+  (N.,VVN.Q.)^  »  (N.,f)^ 
j  11  h  j  11  h  3  h 

+  {G(V'^N.,7N.q.)^+  (N.,V*7N.q.)^  -  (N  .  ,Lq)  }  ,  j  -  1,2 . m  ,  (6) 

jiinj  ixn  j 

where  q^  is  the  value  of  the  asymptotic  solution  at  the  i-th  nodal  point. 
Here  the  terms  included  in  the  first  line  coincide  with  the  SFE  scheme, 
eq.  (5),  while  the  terms  in  the  second  line  (in  curled  brackets)  represent 
the  correction  term  which  is  the  essence  of  the  present  AFE  scheme.  This 
AFE  scheme  can  also  improve  the  pointwise  error  estimate  of  the  SFE  scheme 
by  a  factor  of  the  0(C^''‘^)  [4]. 

Usually  the  finite  element  solution  supplies  values  everywhere  inside 
the  elements  via  the  interpolation  (4) .  Applying  the'  same  interpolation 
using  instead  of  will  not  give  good  results  within  the  elements 

especially  when  there  are  no  nodal  points  inside  the  boundary  layers.  One 
should  use  instead  the  interpolation 

q(x)  =  q(x)  +N^(x) 

which  recovers  the  proper  boundary  layer  behaviour. 

3 .  ASYMPTOTIC  SOLUTION 

The  approximation  q(x,y)  can  be  obtained  by  the  method  of  matched 
asymptotic  expansions  (Cole  [3]) .  Such  approximations  usually  satisfy  the 
boundary  conditions  and  become  increasingly  accurate  as  C  decreases. 

Often  the  error  decreases  like  some  power  of  e  depending  on  the  number  of 
terms  used  in  the  construction  of  the  solution  (see  Table  I) . 

We  first  construct  the  zeroth  order  asymptotic  approximation  from  the 
outer  solution  q  ,  where 


and  from  the  boundary  layer  q^j  satisfying  a  one-dimensional  boundary 
layer  equation  in  the  direction  normal  to  the  boundaries  (where  the  flow 
exits  the  computational  region) .  A  complete  zeroth  order  solution  should 
also  include  corner  regions,  tangential  regions,  and  boundary  layers 
developing  from  discontinuities  in  boundary  conditions.  Restricting  our¬ 
selves  for  now  to  continuous  boundary  conditions  and  constant  velocity 
field  not  tangent  to  any  boundary,  we  find  that  only  the  exit  corner  layer 
has  to  be  included. 

Equation  (7)  implies  that  the  solution  remains  constant  along  stream¬ 
lines,  consequently  it  carries  with  it  the  q-values  entering  the  computa¬ 
tional  region.  The  difference  between  these  values  and  the  values  encount¬ 
ered  at  the  exit  forces  the  boundary  layers. 

We  consider  the  flow  in  the  unit  ^quare  with  vertices  (0,0) ,  (0,1) , 
(1,0) ,  (1,1) .  Let  the  components  of  V  be  u  and  v  (both  positive)  and 

the  differences  of  the  exit  boundary  x  -  1  and  y  ■  1  be  f(y)  and  g(x) 

respectively.  Then  the  structure  of  the  boundary  layers  will  be 

-  u(x-l)/E  v(y-l)/e  u(x-l)/E  v(y-l)/E  . 

■  f(y)e  +g(x)e  -  pe  e  ^  ,  (8) 

where  p  -  f(l)  ■  g(l)  and  f(0)  -  g(0)  0  by  assumption. 

The  third  term  in  (8)  is  the  corner  boundary  layer  and  it  satisfies 
the  differential  equation  exactly  (in  this  particular  case) .  The  approxi¬ 
mation  q  is  obtained  by  adding  q  from  (7)  to  q  from  (8) .  It  is  a 


It  is  a 


unifomly  valid  approximation  in  the  square,  and  is  easily  adapted  to  par¬ 
ticular  cases.  In  our  example  u  «■  v  ■  /2/2  and  the  outer  solution  is 

=  simi(x-y)  ,  (9) 

corresponding  to  the  boundary  conditions 

q(x,0)  =  siniTx;  q(0,y)  »  -siniiy  .  (10) 

On  the  other  hand  we  take 

q(x,l)  =»  simr(x-l)  +  X,  q(l,y)  *  siniT(l-y)+y  ,  (11) 

giving  rise  to  f(y)  *  y  and  g(x)  ■*  x,  which  satisfy  the  requirement  of 
continuity  on  the  boundary. 

4 .  RESULTS 

We  solved  the  differential  equation  (1)  in  the  unit  square  with  the 
boundary  conditions  (10)  and  (11)  on  a  net  with  6,  12  and  24  equal  inter¬ 
vals  in  the  x  and  y  directions.  We  used  centered  three-point  differences 
for  the  first  and  second  derivatives.  As  we  do  not  have  the  exact  solu¬ 
tion  for  this  problem,  we  used  as  reference  a  solution  obtained  with  the 
AFE  scheme  employing  a  uniform  mesh  of  24  x  24  biquadratic  Lagrangian  ele¬ 
ments  (more  details  appear  in  [5]) .  The  range  of  c  reproduced  in  the 
tables  was  C  >  0.01,  £  •  0.02,  e  «  0.05. 

As  a  seminorm  we  used  a  weighted  average  of  the  absolute  value  of 
the  errors.  A  proper  mix  of  interior  and  boundary  layer  points  was  ob¬ 
tained  by  taking  only  the  25  mesh  points  contained  in  a  square  near  the 
corner  (1,1). 

We  observe  that  (Table  Z)  the  error  in  the  asymptotic  approximation 
A(e,h)  decreases  roughly  like  e.  The  error  in  the  regular  (Table  III) 
numerical  solution  N(c,h)  decreases  very  slowly  with  h  and  increases 
as  £  decreases. 

We  note  first  that  the  errors  in  Table  IX  are  in  all  cases  smaller 
than  the  corresponding  errors  in  Tables  I  and  III.  Moreover  the  analysis 
of  [2]  indicates  that  the  error  in  the  improved  solution  B(e,h)  should 
decrease  with  £  and  h  (Table  II)  .  In  fact  it  should  be  proportional 
to  the  product  of  the  previous  errors.  Table  IV  presents  the  ratio 
K(£,h)  »  B (£ ,h) /( {A(e ,h) N(e ,h) ) ) .  We  expect  it  to  approach  a  constant, the 
fact  that  the  values  in  Table  IV  are  not  far  from  unity  maJces  the  booster 
method  quite  attractive. 

We  note  that  the  centered  method  is  quite  useless  by  itself.  Our 
experience  shows  that  the  Booster  Method  worics  equally  well  with  other 
schemes  and  other  discretization  methods. 
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Lectures  on  Spectral  Methods  for  Turbulence  Computations. 

S.  A.  OnszAC  (*) 

Monndrliuarliii  TimlilKlr  of  Trchiiolof/ii  -  Cuiiihriihjr,  MA  02I3U 


1.  -  Introduction  to  spectral  methods. 


Spectral  methods  are  based  on  representing  the  solution  to  a  problem  as 
a  truncated  series  of  smooth  functions  of  the  independent  variables.  Whereas 
finite-element  methods  are  based  on  expansions  in  local  basis  functions,  spectral 
methods  are  based  on  expansions  in  global  functions.  Spectral  methods  are 
the  extension  of  the  standard  technique  of  separation  of  variables  to  the  solution 
of  arbitrarily  complicated  problems. 

Let  us  begin  by  illustrating  spectral  methods  for  the  simple  one-dimensional 
heat  equation.  Consider  the  mixed  initial-boundary- value  problem 


(1.1) 


t)  _  jj-  t) 

et  ■■  ex‘ 


(0<  ar<  ,x,  t  >  0), 


(1.2)  t»{0,t)  =«(.x,/)  =  0  (f>0), 

(1.3)  u(x,0)=f(T)  (0<x<.i). 


The  solution  to  this  problem  is 


(1.4)  u(x,  t)  ~  ^  a,(t}  siu  nx  , 

n"l 

(1.5)  a,(f)  = /,exp[— , 
where 


(l.«) 


-El 


/(j)  sin  7tJ-dr 


arc  the  coeflicients  of  tlie  Fourier  sign  series  c.\pnIl.^ioll  of  /(x). 
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A  sjicctrnl  :i]i|iro.\jiii:itiiiii  t<i  ( 1 .1 )-( 1  is  "(ittiMi  t)\  siiiiiily  Iruiicatiii^'  ^ 

(1.1)  tr.  I 

A  ' 

(l.T)  «v(.i',  t)  —  2  f(„(0  sin  )i.f 

H  -I  < 


iHKl  it‘i)l:uin}:  (l-.j)  by  the  evolution  eqiiiition 

(1.8)  =  —  A'»*o,  (u  =1, ....  *V) 

Avitli  the  initial  condition.s  <»„(0)  =  /„  («  =  l, ....  A'). 

The  spectral  approximation  (1.7),  (1.8)  to  {1.1)-(1.3)  is  an  exceedinply  "ood 
approximation  for  any  time  t  greater  than  zero  as  N-*  oo.  In  fact,  the  error 
u{x,  t)  —  M.v(j‘,  t]  satisfies 

A 

(1.9)  tt(x,  t)  —  «.v(J',  t)  =  2  .fn  pxp  [— A'w*t]  sin  wj- =  0(exp  [— TlA'*/]) 


for  any  <>  0.  In  contrast  to  (1.9),  finite-difference  approximations  to  the  heat 
equation  using  N  grid  points  in  x  lead  to  errors  that  decay  only  algebraically 
with  A  as  A  oo.  Furthermore,  this  spectral  method  for  the  solution  of  the 
heat  equation  is  eflftciently  implementable  by  the  fast  Fourier  transform  (FFT) 
in  O(AlogAr)  operations. 

There  are  several  significant  difficulties  in  extending  the  simple  spectral 
method  employed  for  (1.1)-(1.3)  to  more  general  problems.  Among  these  dif¬ 
ficulties  are  those  caused  by  imposition  of  nontrivial  boundary  conditions, 
nonlinear  and  nonconstant  coefficient  terms,  and  complex  geometries.  These 
difficulties  and  their  solutions  will  be  discussed  below  (see  al.so  [1,2]). 

The  Fourier  series  (1.4)  converges  fast  if  tt(x,  t)  is  infinitely  differentiable 
and  u(t,  t)  satisfies  the  boundary  conditions 


;  (i.u)) 

> 


c’"m(x,  t) 
cx*» 


=  0 


(X  =  0,  .t) 


for  all  nonnegative  integer.^  n.  Under  these  conditions,  the  error  after  A'  terms 


I 


f.v(x,  f)  ==  tt(.r,  0  —  2  ®"(0  ><i>'  »  r 


goes  to  Zero  uniformly  in  .r  faster  than  any  jxower  of  1;A'  as  A' — oo.  On  the 
other  hand,  if  «(x,  t)  is  not  iitfinitely  difTeientiahle  or  if  any  *if  the  coiiditinns 
(1. 10)  is  violatx'd,  then  ry(.r, /}  =  0(1. '.V'')  ns  A’  -  oo  for  sonn*  j).  For 


i.i;i  iiiir.s  i>\  rii'r.i'i  It  M.  mi.iih>ii»  iok  i  i  itiii  i.k.m  i;  <'ii\im  m  iiivt 
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fXMlIlplf, 

(1.11) 

l)Ut  ttio  erroi-  iiicurwl  l>y  trunc-iitiiij:  aftiT  X  tonus  is  of  older  IjX  for  any 
fixed  j  ,0  <  X  <  :i.  Furtliennore,  the  eoiiver^tenee  of  (1.11)  is  not  iiiiifoi  ui 
in  .r;  (1.11)  exhibit.s  Gibb.s’  jdienonieuon,  iiiimely 

f  v(|/*V)  =  0(1)  (*v  oo,  I  fixed)  . 

For  any  fixed  X,  there  are  points  x  at  which  the  error  after  X  terms  of  (1.11) 
is  not  small.  The  poor  convergence  of  (1.11)  is  due  to  the  x-iolation  of  (1.10) 
for  n  =  0. 

More  generally,  most  eigenfunction  expansions  of  a  function  /(x)  converge 
faster  than  algebraically  (i.e.  the  error  incurred  by  truncating  after  X  terms 
goes  to  zero  faster  than  any  finite  power  of  l/.y  as  X-*  oo)  only  if  /(x)  is  infi¬ 
nitely  differentiable  and  /(x)  satisfies  an  infinite  number  of  special  boundary 
conditions.  For  example,  the  Fourier-Eessel  expansion 

f(Jf)  =  ^a,J,{XnX)  (0<.r<l), 


I(-l)' 


sin  C-’ii  -•  1  ).<■ 
'Jn  -j-  I 


(0 


■T)  , 


where  Xn  is  the  «-th  smallest  root  of  J,{X)  =  0,  converges  faster  than  algebrai¬ 
cally  only  if  /  is  infinitely  differentiable  and 

(1.12)  =  1 

for  k  =  0, 1,  2, .... 

Wlien  a  spectral  expansion  converges  only  algebraically  fast,  spectral 
methods  based  on  these  eigenfunction  expansions  cannot  offer  significant 
advantages  over  more  conventional  (finite-difference,  finite-element)  methods. 
Eigenfunction  expansions  of  this  kind  should  not  normally  be  used  unless 
tlie  boundary  conditions  of  the  problem  imply  all  the  extra  boundary  con¬ 
straints  like  (1.10)  or  (1.12).  For  example,  if  periodic  boundary  conditions 
are  compatible  with  the  differential  equation  to  be  solved,  complex  Fourier 
series  are  suitable  to  develop  efficient  spectral  approximations. 

In  the  development  of  spectral  methods  for  general  problems,  it  is  important 
tliat  the  rate  of  convergence  of  the  eigenfunction  expansion  being  used  docs 
not  depend  on  special  properties  of  the  eigenfunctions,  like  boundary  conditions, 
but  ratlier  depend  only  on  the  smoothness  of  the  function  being  expanded.  Of 
course,  if  the  solution  to  the  problem  being  solved  is  not  smooth,  one  should 
not  expect  errors  that  decrease  faster  than  algebraically  with  IjX  when  global 


'-.Vi 


f 
¥ 


K 


Vi 


no 


S.  A.  l>liS7vi: 


(■\])aii8i<>iis  art*  used.  Kastor  than  alfri'biaii-  i  atcs  (if  ((iiiA  cn’^ciicc 
may  be  acbicvotl  for  tboso  jirobk-ms  by  citlior  iiatchiiiK  ll;c  solution  at  discon¬ 
tinuities  or  lire-  and  iiosf-jiroeessiii}:  tf  tlie  solution  (ste  [Itj). 

Then-  is  an  easy  way  to  ensure  that  the  rate  of  conver'reme  of  a  speetral 
expansion  of  a  funetion  f(x)  depi-nds  only  on  the  smoothness  of  /(.r),  not  on  its 
boundary  properties.  The  idea  is  to  expand  in  terms  of  suitable  elusi-es  of 
orthogonal  jiolynomials,  including  Chebyshev  and  Legendre  polyncinials  for 
all  those  problems  iji  which  constraints  like  (1.10)  ajid  (1.12)  are  unrealistic. 
These  polynomial  expansions  avoid  all  diificulties  associated  with  the  Gibbs 
phenomenon  provided  the  solution  f{T)  is  .emcoth. 

From  the  mathematical  point  of  view,  the  classical  orthogonal  polynomials 
are  eigenfunctions  of  singular  Sturm-Liouville  problems.  It  is  not  hard  to 
show  [1]  that  expansions  using  eigenfunctions  of  such  singular  Sturm-Liouville 
problems  converge  at  a  rate  that  depends  only  on  the  smoothness  of  f(x),  in 
contrast  to  eigenfunction  expansions  based  on  nonsingular  Stiurm-Liouville 
problems  that  lead  to  additional  boundary  constraints  like  (1.10)  on  f(x). 

These  results  for  orthogonal  polj'nomial  expansions  are  easily  demonstrated 
in  the  case  of  Chebyshev  polynomial  expansions.  The  n-th-degree  Chebyshev 
polynomial  T„(x)  is  defined  by 


(1.13) 

Therefore,  if 

(1.14) 
then 

(1.15) 


T,(cos  6)  =  cos  nQ  , 


f(x)  =  2  a.3’.(®) , 


—  /(cos  0)  =  2  cos  . 


Thus  the  Chebyshev  polynomial  expansion  coefficients  a,  of  /(r)  are  just  the 
Fourier  cosine  expansion  coefficients  of  the  even,  periodic  function  g(0).  A  simple 
integration-by-parts  argument  then  shows  that 


n*a,->  0 


(»  ->oo) , 


provided  gifi)  (or,  equivalently,  /(x))  has  p  continuous  derivatives.  Since 

[/(^•)  -  fl»r-(j?)j<  Z  [“"I  (l-*'l<i)» 

it  follows  that  the  rate  of  convergence  of  (1.14)  is  faster  than  algebraic  if  / 
is  smooth. 
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Ill  sumiiiiiry,  spoi-tral  (“Xjiaiisions  .slimiltl  lie  niado  u.siii;:  scrii's  of  ortli(i{;oiial 
])olynomials  unless  tlu-  hoiindary  conditions  of  tlic  prolilcin  arc  fully  conipatildc 
with  some  other  class  of  ei;fenfuiic!ioiis.  In  iiraclice,  Clieliyslicv  and  Legendre 
polynomial  t'xjiatisions  tire  reconi mended  for  most  a]>]ilieations,  supplemented 
by  Fourier  series  and  surface  harmonic  series  when  hoiindary  conditions  iierniit. 

Another  difficulty  with  general  kinds  of  spectral  metheds  is  their  a]iiilication 
to  problems  with  nonlinear  and  nonconstajit  t  oeflicient  terms.  Before  explaining 
the  solution  to  this  problem,  let  us  illustrate  the  difficulty. 

Suppose  we  wish  to  solve  the  partial  differential  equation 

(1.16)  ~  =  jr(u,  u)  +  seu , 

C  *• 

where  u  =  «(jc,  <)  and  is  a  bilinear  (nonlinear)  operator  that  involves  only 
spatial  derivatives  and  .Sf  is  a  linear  operator  that  involves  only  spatial  deriva¬ 
tives,  The  operators  J/"  and  may  depend  on  both  *  and  t.  A  spectral  method 
for  the  solution  of  (1.16)  is  obtained  by  seeking  the  solution  as  a  finite  spectral 
expansion : 

(1.17)  «(x,  <)  =  2  <».(0V.(*) » 

•  “1 


where  we  assume  for  now  that  y',(x)  (1  <n  <  oo)  are  a  complete  set  of  orthogonal 
functions.  If  we  introduce  the  re-expansion  coefficients  e,.,  and  so  that 


*1-1 


and  equate  coefficients  of  Vs,(*)  (n  =  1,  in  obtain 

da  ^  ^ 

(1,18)  -77  =  2  2 (n  =  l,...,A). 

n*  m*!  fH— I 

Equations  (1.18)  are  the  spectral  evolution  equations  for  the  solution  of 
(1,16).  They  have  one  very  serious  drawback.  In  general  c,*,  and  d,«  are 
nonzero  for  typical  n,m,p,  so  that  evaluation  of  da/dl  from  (1.18)  for  all 
n=  1,...,A''  requires  0(A'*)  arithmetic  operations  for  the  bilinear  term  and 
O(N^)  operations  for  the  linear  term.  Thus  solution  of  (1.18)  requires  order 
operations  per  time  step.  Since  operational  siiectral  calculations  now  involve 
A''>10*,  the  computational  cost  of  the  direct  solution  of  (1.18)  is  prohibitive 
(even  if  only  linear  terms  arc  present). 


si 


Cl 

I 


i 


r^rwruTv  )L«  K.*  )C".v.v  ".--w 
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Ill;'/,  vi; 


Till-  luoliloin  li«T<‘  is  oiu-  Ilf  <'<iin|Mitiitii>iiiiI  cinniilfxil v.  l''iiiitc-iliili'H’iiri' 
iml hulls  fur  till'  siilittiiiii  iif  (l.lf!)  on  X  ^rriil  jioiiits  mny  ri‘nnir<‘  only  oiilor 
X  o|MT;if ion.s  ]K‘r  tiini-  .sfojt.  If  llic  s])citi’:il  inofhoil  ri'ally  ri‘i|uiirs  onh  r  .V’ 
oiionitioiis  j)cr  tiiiu‘  sti'j).  it  cannot  l•oInJ)cl<‘  Avhcn  -Y  is  lar^'c. 

Anotlicr  cxainjilc  illustrating;  the  <-oinpntatiiinal  complexify  of  sjM-itral 
nictlioils  is  ;rivcn  by  the  nonlincc.r  clilTiision  cipiation 


(I.IH) 


ct 


C*H 

=  oxp  [((.]  (x,t). 


If  Avc  seek  the  .solution  a.s 


(1.20) 


t<(x,  t)  =  "2,  o„{t)ip.(x) 

n—l 


in  terms  of  the  orthonoi*mal  functions  y.\{x),  then 
da  r  r  ^  1  ^ 

(1.21)  “df  J  I  ^  <»m(t)V«(j:)]  2  «pV’I(J’)d.r 

for  n— These  evolution  equations  for  {a,(f)}  have  an  exponential 
degree  of  computational  complexity  as  they  are  expressed  as  an  integral  func¬ 
tional  of  {a,(l)}. 

The  solution  to  the  problem  of  computational  complexity  is  to  use  the 
author’s  transform  methods.  Let  us  illustrate  the  technique  for  a  pseudo- 
spectraJ  (or  collocation)  approximation  to  (1.19).  First,  we  introduce  X  suitable 
collocation  points  j,,  jr,, ....  afy  lying  within  the  computational  domain.  Then 
the  approximate  solution  (1.20)  is  forced  to  satisfy  the  partial  differential 
equation  (1.19)  (or  its  boundary  conditions)  exactly  at  these  discrete  points 
at  every  time  1.  More  specifically,  the  following  three  steps  are  done  at  each 
time  step  (: 

i)  Determine  A*  coefficients  a„(t)  (n  —  1, X)  so  that 


(1.22)  u{x„  f)  =  2  n«(nvi,(x,) 
ii)  Evaluate  u,Jx,,t)  by 

(1.23)  ««(J’,.  1)  =  2 


ij  =  1 . A)  . 


(1  =  1 . -V). 


(1.21) 


iii)  Finally,  evaluate  cM(x,.f),cf  by 
c«(x,,  f) 


ct 


cxp[u(.r,,  t)]i/„(,r,,  f) 


(/  =  1,....  -V) 


and  Tiiarch  forward  to  the  next  time  .step. 


I.liKKKS  ii\  >ci;<ri( \l.  UKIiliiKS  foti  l(  ttlt(  I.KNCI.  CoMITTA  IION  < 


11.3 


'I’lic  i<l(M  of  fli<‘  |)si‘iiilo.s{M‘i-traI  (raii.sforin  nu'tlKxl  lx*  ix’staled  as 
3'ran>forni  frcoly  hotwooii  physical  (.r,)  ainl  spectral  (<i,)  representations,  evaluat- 
in;:  each  term  in  whatever  representation  liiat  term  is  most  accurately,  aiul 
simply,  evaluated.  Tims,  in  (1.21),  we  evaluate  exp  [u]  in  the  physical  represen¬ 
tation  while  we  compjite  in  the  spectral  repre.sentation  by  (1.22)  because  it 
is  most  accurately  done  there. 

It  should  be  apparent  to  the  reader  that  pseudospectral  transforni  methods 
can  be  applied  to  any  problem  that  can  be  treated  by  finite-difference  methods 
rejiardless  of  the  technical  complexity  of  nonlinear  and  nonconstant  coefficient 

terms. 

For  the  expressions  of  interest,  eomputation  of  derivatives  of  a  A’-term 
si)ectral  expansion  require.s  order  A'  arithmetical  operations.  For  the  Fourier 
series  (1.7),  this  fact  is  obvious; 


d  '  ” 

ojr  ,_i 


d*  * 

;»  »  2  nx  =  —^  rt*a,  sin  nx  . 

U-T* 


For  the  Chebyshev  polynomial  expansion  (1.14),  the  computational  complexit.v 
of  differentiation  is  a  little  less  apparent.  Since  T«(cos  6)  ==  cos  nd, 


n-f  1  n  — 1  c,  ' 


(n>0), 


where  c,  =  2,  c»  =  1  (n>l)  and  T'  =  Tl,  =  0.  Therefore,  if 

d 


.  2  «.  r,(x)  =  2  > 

u-r 


then 


y  r  T  T'  T'  1 

-*  I  o,  t:(x)  =  V  r.h,  =  2  • 


Equating  coefficients  of  r’(x)  for  n  =  l,...,A’'  -f  1  gives  the  recurrence  re¬ 
lation 


(1. 


<',-1^-1  - 
*  =0 


(l<»i<A'), 
(»i>A’) . 


The  solution  of  (1 .25)  for  h,  given  <i,  requires  only  order  A’  arithmetic  operation-. 
•Similar  recurrence  relations  can  be  obtained  for  diffeientiation  of  spcctr.il 
scries  based  on  other  sets  of  orthogonal  polynomials  and  functions. 
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Ill  tlic  <;iso  of  Fouvior  sorios,  llio  traiisfuriii  (l.T)  ami  itN  inviT.'O  ran  1m! 
i‘i»i)]iuti‘il  ill  0(.V  lov'j -V)  operations  if  X  - tlie  fa>t  Foiiiii'i'  traiis- 
forni.  However,  most  of  llie  eompiitatioiial  eflicieiiey  of  transform  metliinls 
comes  not  from  the  FFT  but  from  the  separability  of  multitlimensional  trans¬ 
forms.  Thus  a  thnH‘-(liinensioiuil  discrete  Fourier  transform  eaii  he  expri'ssed 
as  three  oiie-diineiisiomil  Fourier  transforms 


1 1 1"'''  “• "  ""  H  (t  +  r + t)]  = 


=  y  exp  ['J.ri/m  'J]  ^  exp  [S.tiA  w  'A'J  ^  a(j,  I,  I)  exp  [•J:tilplL] . 

J-»  l-»  1-0 


The  left-liand  side  of  (1.26)  requires  roughly  {JKL)*  operations  to  evaluate 
at  all  the  points  0<»t<t7,  0<n<  A',  0<p<i,  On  the  other  hand,  even 
without  the  FFT,  the  right-hand  side  of  (1.26)  requires  only  about  (JKL)- 
'(J  -r  K  —  L)  operations  to  evaluate  at  all  the  points.  When  the  FFT  is  applied 
to  the  one-dimensional  transforms  on  the  right-hand  side  of  (1.26),  the  number 
of  operations  necessary  to  evaluate  (1.26)  is  reduced  further  to  (JKL)- 
•(log,./ 4-  log, A 4-  logtZ)  if  J,  K,  L  are  powers  of  2. 

Spectral  approximations  to  general  boundary-value  problems  lead  to  full 
2fxy  matrix  equations  for  the  N  expansion  coefficients  a,.  It  would  seem 
that  solution  of  these  equations  requires  OiK*)  arithmetic  operations,  while 
storage  of  the  matrix  requires  0(A*)  memory  locations.  Since  tj'pical  problems 
now  involve  A '*<10*,  the  direct  solution  (or  even  the  direct  formulation)  of 
such  problems  would  seem  unworkable  now. 

Consider  the  solution  of  a  general  linear  differential  equation  Lu  =  f.  Let 
a  A-term  spectral  approximation  to  this  problem  be  given  by 

(1-2:) 

where  /.>•  is  a  suitable  A-term  approximation  to  /.  As  mentioned  several  times 
earlier,  the  matrix  representation  of  (1.27)  is  generally  a  full  Ax  A'  matrix,  so 
that  direct  solution  of  (1.27)  by  Gauss  elimination  methods  would  require 
order  storage  (for  the  matrix  representation  of  Z,^)  and  order  A*’  arithmetic 
operations. 

Here  we  shall  describe  a  method  that  permits  the  solution  of  (1.27)  using 
order  A'  storage  locations  with  the  number  of  arithmetic  operations  of  order 
the  larger  of  A'  log  A'  and  the  number  of  operations  required  to  solve  /,«  =  f 
by  a  iir»t-ordcr  fuiite-diflerencc  method.  The  important  conclusion  is  that 
xpfctral  mflhodH  for  grneral  problems  in  grnrrnJ  gromrtries  ran  be  implrmeytied 
efficienthi  tcilh  operation  costs  and  storage  not  much  larger  than  that  of  the  simplest 
finite-difference  approximation  to  the  problem  irith  the  same  number  of  degrees  of 
freedom.  Since  s])cctiiil  methods  rcquiri-  mciiy  fewer  dcgi'c*<s  cf  frccdc'm  to 


^  ^ 
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”i  \  ('ll  iicfuriicv  (or,  li«';irl\  «-(|lli\  ;i  lout  1  \ .  n|hm  (  r:il  mcl  !;<  «!>  ;  i  liic\  «•  1 11 II i  ll 
lii;:ln'i  iiccuri'.cy  lor  a  iiiiiiil:('r  of  i1<';:i'('cn  of  fiti'il'.  in)  tli!;ii  i'(<niir(il  liy 

liiiilc-oritor  fiiiilf-iliri'croiiu*  api'ro.viin:  tioii.'.  inipoi  li'iil  <  <;ni[iiit;:(ioiial  <111- 
<  i<*iiii<‘.s  ifsult  from  tin-  lu  w  m»-tl;«)<l. 

Tli<‘  i<lca  <if  tli»“  iloratioii  m<-tho<l  is  as  follows;  Suiiiio.-c  we  arc  able  to  coii- 
stiiict  an  aiii>ioxim:!lion  to  flu*  sjioifial  oiu'rator  tliat  lias  tl:<*  following* 
jirojicrtios: 

i)  lias  a  sparse  matrix  representation  so  tliat  it  l  aii  be  represented 
usiii};  only  OfA’)  storage  locations. 

ii)  is  eftieiently  invertible  in  the  sense  that  the  equation 

4p«.v=/.V 

is  sahTible  as  effleiently  as  a  first-order  finite-differenee  approximation  to  the 
problem. 

iii)  approximates  in  the  sense  that 

(1.29)  0<m<]:i-'X.,:;<J/<  oo 

for  suitable  constants  m,  if  as  AT  -*■  oo.  Roughly  speaking,  (1.29)  requires 
that  the  eigenvalues  of  X^pX,,  be  bounded  from  above  and  below  as  A’  -»■  oo. 

We  propose  to  construct  L„  from  X„  by  changing  the  discretization  operator 
either  in  addition  to  or  in  place  of  approximating  the  differential  operator. 
Thus  we  construct  X,„  by  a  suitable  loic-order  fiuHe-differcncc  approximation  to  L. 
A  simple  example  is  given  by  the  second-order  differential  equation 

(1.30)  i«  =  /(r)u'(x)  -f  g(x)u'[x)  -f  h{x)u(x)  =  r(.r)  (0<.r<  2.t) 

with  periodic  boundary  conditions  «(x  -f-  2.t)  =  «(x)  and  /(x)  >  0.  A  spectral 
upproximution  is  approximately  sought  as  the  finite  Fourier  series 

(1.31)  «(x)  =  a^rxp[iLr]  . 

:ti<M 

If  the  Fourier  eceffieients  of  f{x),  g{x),  Ii(x).  r(x)  are  denoted  fo  .  r^, 

respectively,  then  the  spectral  (Galeikin)  equations  for  aie 

(1 .32)  X.„n  =  2  [-  -t-  T>sr,_,  -f  /q  .]«,  =  . 


(Mrurly,  tlio.so  t^quatiiMis  in  a  full  matrix  r^'prosfutaticn  that 

iTfiiiirrs  0{J\-)  storajrr  hniiti<ns  and  0(7i-)  (i|U'i at iiiis  tu  inxi-il. 


lift  >.  V.  MK-7.\‘. 

A  suitaMf  :(pt)roKiin:ito  oin'rator  7-^,,  is  *'(Mi>tni('l<Ml  luin;:  the  colliu-af i(Hi 

jiuints  .r^ -- •Jrrj/A’’  (j  ^  o,  .1 . -V— I),  \\li(*r<‘  ‘2J\\  In  tlu*  plivsical  .x|»ar<* 

I'cpicsontation,  wo  uso  tlio  fiijito-dilloriMuo  approxiiiiat ion 

(l.;wi  X„1<'„=:  /(.r,)  ^  ^  ’ 

avIh'io  Uj=  h(j',)  and  Aj"  =  2.^/A■.  Obviously,  X,„  is  sjjarso  and  ifHcieiitly  in¬ 
vertible.  To  verify  (1.29)  ue  u.se  the  following:  elementary  argument  (that 
may  be  made  more  rigorou.s  but  no  more  correet  by  more  involved  WKll-like 
arguments).  If  A  is  an  eigenvalue  of  then  there  exists  .a  function  u{.r) 

gueh  that 

(1.31)  i.,«  =  . 

If  u(x)  is  a  smooth  function  of  x  (in  the  limit  A’-vcx)),  then  both  and 
i.jM  should  be  good  approximations  to  Z«(x),  so  (1.34)  implies  /.~1.  On  the 
other  hand,  if  «(.r)  i.s  a  highl}’  oscillatory  function  of  .r  (in  the  limit  X  ^  oc), 
then 


(.V  oc)  . 


Therefore, 


X.,  «  —  /(•>•») 


—  2«, -f  M,_, 


and,  if  transform  (pseudospectral)  methods  are  used  to  evaluate  L,^v, 


so  (6.18)  gives 


(1.38)  /(.r,)  2  (— l’)«tCxp[tl-r,]~/./(.r,) 

U-.<r  l-a./  ) 


The  eigenfunction.s  of  (1.38)  are 


«>  =  <-xp[/7.r,] 


iW.  'in 


and  the  associated  eigenvalue  is 


(ly  Art' 

1  sin*  1 7 \f 


-'.v.v.v.v.' •, 
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SiiKc  n'  <  A'  willi  ^ -V  =  .T  A.i .  Mc  (ililiiiii 


'I'lais  (1. liiilds  with  w  =  l  iiiul  J/ =  .1-;  I  '5:2."). 

Tlicrc  arc  scvcial  cxt«'iisi<iiis  of  the  above  iiietlKiil  for  eniih;trtietii!e 
tliat  are  iinpoi  tant  in  ]traeti<e.  First,  in  the  ease  of  Clu'bysliev  .spef  tral  niethmls, 
if  is  atipropriate  to  «onstruet  iisiii}:  tinite-dillereiue  approximations  based 
on  tlie  collocation  points  j’,  =  cos, -rj/A’.  In  this  ease,  the  operator  bounds 
(I.2t>)  continue  to  hold  with  21—  ‘2.5,  m  =  1  for  a  wide  variety  of  operators  L. 
Second,  higher-order  equations  are  best  treated  by  writing  them  as  a  system  of 
lower-order  equations.  Thus  direct,  consti'iiction  of  for  i  =  V‘  gives 

Howevei',  if  we  introduce  t  —  V-u  and  define  the  second-order  operator  K  by 


then  direct  construction  of  A'„  as  a  finite-difference  operator  gives 


Third,  odd-order  operators,  initial-value  problems  and  problems  of  mixed  type 
are  best  treated  by  constructing  on  a  grid  that  is  roughly  50%  finer  than 
that  used  in  construction  of  by  collocation.  In  this  case  the  spectral  bounds 
(1.29)  with  J/<2.5  continue  to  hold  for  most  problems.  For  example,  the 
operator  c/eJa;  with  periodic  boundary  conditions  has  spectrum  t'A,  while  its 
centered  finite-difference  approximation  has  spectrum  t  sin  {kAx)lAr,  so 

'  =  0(t  Ax/sin  I- Ax) , 


which  is  unbounded  for  |i  AjI<:r,  but  bounded  by  4.i/3V3  ^  2.4  if  |A  Axl<  2.t,'3. 


2.  -  .Applications, 


2'1.  Introduction.  -  Over  the  last  few  years,  there  has  been  progress  in 
understanding  fundamental  nonlinear  jirocesses  in  she.ar  flows.  In  this  section, 
I  shall  survey  some  result.s  that  have  emerged  from  numerical  studies  of  tran- 
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.silidii  and  I  mludciicv.  i  sliall  r«-vii'\v  for  you  tlncf  dillVi  ciil  a.''pfil>  «il'  IIk'nc 
](r(d)lcins.  First,  I  sliall  suiniiiari/f  rl•s^llts  on  tin-  basic  i iivla bilil ics  that  scciii 
to  be  icsiimisiblc  for  till*  onset  of  ebaos  in  these  Hows.  These  instabilities  apiiear 
to  be  ntliver.sal  in  eliaracter  and  may  e\)>lcin  many  of  the  nnifyitij:  features  of 
transition.  Second.  1  Khali  {five  koiik'  e.vam]iles  of  )iro;.TesK  in  the  numerical 
simulati'.m  of  bi;jh-l{eyiwdds-niimber  Hows.  Finally,  I  will  cive  a  synojisis  of 
new  idiMs  for  .siib^rid  .scale  c!().snre.s  of  buce-Iievnolds-munber  tiirbaleme. 

Full  details  of  the  ideas  discussed  here  tire  civen  in  the  references. 

2'2.  .1  IranitUional  in»tability.  -  The  processes  by  which  laminar  flows  under¬ 
go  transition  to  turbuleuee  remain  basically  unsolved.  However,  recent  numeri¬ 
cal  studies  have  provided  some  insights  into  transition,  including: 

2‘g.l.  Xonclassical  character  of  transitional  instabilities. 
The  primary  linear  (exponential)  instability  of  classical  plane  parallel  shear 
flows  Avith  noninflectional  velocity  profiles,  as  described  by  the  Orr-Somnierfeld 
(or  related)  equations,  is  much  too  weak  to  describe  transition.  For  ex¬ 
ample,  linear  instability  of  plane  Poiseuille  flow  [V(z)  =  l  —  z^,  1*|  <  l)  oc¬ 
curs  for  Reynolds  numbers  .R,,>5T78,  while  Squire’s  theorem  implies  that 
the  critical  disturbance  is  two-dimensional.  The  fact  that  this  instability  is 
induced  by  a  subtle  interplay  of  viscosity  and  shear  implies  that  its  growth 
rates  are  quite  small  on  convective  time  scales.  For  example,  the  most  rapidly 
growing  exponential  mode  of  the  Orr-Sommerfeld  equation  is  obtained  at 
—  48000;  its  growth  rate  is  only  0.0076;  it  is  so  feeble  that  perturbations 
grow  by  a  factor  10  in  a  time  of  about  300,  in  which  time  a  point  on  the  center- 
line  moves  about  150  channel  widths.  In  contrast,  transition  is  observed  to 
occur  explosively  over  a  few  channel  widths  at  Rejmolds  numbers  as  low  as 
roughly  1000.  A  transitional  instability  that  affects  noninflectional  plane 
parallel  shear  flows  must  have  a  characteristic  convective  time  scale. 

2’2.2.  Three  dimensionality  of  transition.  Two-dimensional  fluids 
do  not  appear  to  exhibit  the  kind  of  strong  chaos  that  is  characteri.stic  of  tur¬ 
bulent  shear  flows.  In  thermal  convection,  Curby  et  al.  [3]  show  that  two- 
dimensional  flows  do  not  appear  to  act  in  a  strongly  chaotic  way,  but  three- 
dimensional  flows  may  be  strongly  chaotic  at  large  enough  Reynolds  number. 
Even  for  inflectional  free  shear  flows,  in  which  there  are  strong  inviscid  tw  o- 
dimensional  instabilities,  Brachet  and  Obszag  [4]  show  that  the  flows  that 
dcA'clop  from  two-dimensional  finite-amplitude  disturbances  are  not  strongly 
chaotic,  in  contrast  to  the  flows  that  develop  three-dimensionally. 

2‘2.3.  Instability  of  two-dimensional  nonlinear  travelling 
waA'es.  Perhaps  the  simplest  instability  that  has  the  character  of  a  transi¬ 
tional  instability  i.s  the  linear  three-dimensional  instability  of  two-dimcnsioiuil 
finite-amplitude  flows.  Orszag  and  Kells  (5]  and  Oilszag  and  P.vteu.v  |(>J 


Fig.  1.  -  Streamline*  o(  the  steady  (stable)  finite-amplitude  two-dimensional  travel¬ 
ling  wave  for  plane  Poiseuillc  flow  at  Jf  =  40(Ki,  plotted  in  the  rest  frame  of  the  wave 
(from  [6]). 

show  how  such  an  instability  fits  the  basic  features  of  transition  in  classical 
shear  flows,  including  their  convective  grov  th  rates,  inherent  three-dimension¬ 
ality,  onset  at  Reynolds  numbers  in  accord  with  experimental  obserA'ations 
and  flow  features  in  accord  with  early  transitional  flows.  Tliese  instabilities 
have  been  analyzed  both  by  direct  numerical  simulation  of  the  CA-olving  tliree- 
dinien.sional  flow  and  by  a  linear  perturbation  analysis  cf  the  nonparaliel  two- 
dimensional  (nonlinear  travelling  wave)  flou.  In  lig.  1,  we  .sliow  the  stream¬ 
lines  of  a  typical  two-dimensional  ba.se  state  (here  for  plane  I’oiseuille  flow  at 
Jl  —  4000).  The  nonparallel  character  of  the  base  flow  leads  to  considerable 
complication  in  its  linear  stability  analysis  (see  [6]  for  the  formulation  of  these 
large-matrix  eigenvalue  problems).  A  topic  of  much  current  research  interest 
i.s  the  development  of  efficient  numerical  methods  for  finding  eigenvalues  of 
the  very  large  matrices  encountered  in  problems  of  this  sort.  In  fig.  we 
give  a  stability  di.agram  for  tlii.s  transitional  instability;  here  we  i>lof  contours 
of  constant  growth  rate  as  a  function  of  the  amplitude  of  the  two-dimensional 
base  state  and  the  Reymilds  number.  The  growth  rates  of  tin's  instability  are 
orders  of  magnitude  larger  than  those  of  (Irr-isommei feld  modes.  The 
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Fig.  2.  -  Contours  of  constant  growth  rate  (labelled  by  growth  rate)  as  a  function 
of  It  and  the  amplitude  of  the  background  two-dimensional  nonlinear  wave  (see  right- 
hand  scale). 

development  of  this  three-dimensional  secondary  instability  seems  to  be  con¬ 
sistent  with  available  experimental  data  on  early  transitional  flows.  In  fig.  3, 
we  compare  contours  of  the  x  velocity  at  the  so-called  one-spite  stage  of  tran¬ 
sition  in  plane  Poiseuille  flow  obtained  a)  experimentally  by  Nishioka,  Iida 
and  Ka>'Batashi  [7]  and  b)  numerically  by  Kleisee  and  Schitjax.v  [8].  The 
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Fig.  3.  -  Contoiir.s  of  x  velocity  in  the  (x,  y)-planc  at  the  one-spike  stage  in  the  labora¬ 
tory  experiments  of  Xishioka  et  al.  (a))  and  in  the  numerical  simulation  of  Kleiserand 
Schumann  (6))  (from  [8]). 


I.l  !•  1 1  iM  '  '  i\  'iTi  1 1:  M.  Ml.  I  liiiii-  lilt:  n  1.1  M  I  '  >  I'l  ii  i  1 1  !■  'N  ' 


121 


li(i\v.>  tliiit  ili'M'liij)  ti'iiiu  till'  initial  liiirar  iii>laliil>t  \  a|i|ii'.  r  tu  liail  ilii'iitli  in 
rlian.s  ami  ( inliuli'ln-i-  amt  imt  lit  .salliralf  ilitn  nrili'iril,  laminar  (Inu  vl.ili 
Similar  in>taliilit ir.'N  liavi-  lll•l‘ll  iniiml  in  liiiiimlai>  lavri',  iilain-  (  nm  iic  ilnu, 
liijic  ^ni^<‘uill<•  Ilow  (sit  [t)|)  ami  in  fn  n  slii'i  r  tlnw  s  (m c  ]  I  j). 


2'2.l.  Coin  iii‘t  i  t  ion  botwoiMi  t\\  o-iliini'iisional  iiaiiiii;:  and  tlin-c- 
ilinieii.siona  1  i  nsl  a  hili  t  it-.s.  IntliTtioiiiiI  frci*  shear  llo\\  s.  like  ini  sin  n  la  \  ci  s 
uml  jet.s,  are  inviseiiily  unstable  to  two-ilimeiisiniial  ilistiirhanees.  Snuiri  ‘s 
thoorein  iniiilie.s  that  these  in.stabilities  are  stron;;est  when  tivn-diiiiensionai; 
when  these  two-dimensional  instabilities  evolve  in  time,  they  saturate  into 
ordered  laminar-tlow  states  eharacterized  by  lartre-seale  vortieal  How  structures. 
These  vortieal  flows  may  themselves  be  unstable  to  subharnionii  (pairinj;) 
in, stabilities,  in  which  two  (or  more)  vortices  are  ])aired  and  j;enerato  a  new 
larper-scale  vortex  motion  (9].  In  these  flows,  the  tliree-dimensional  instability 
discussed  above  is  also  present  [10],  but  it  is  not  necessarily  stronger  than  the 
pairing  instability.  However,  the  three-dimensional  secondary  instability  i.s 
effective  at  much  smaller  spanwi.se  spatial  scales  than  is  the  inviscid  primary 
instability  and  seems  to  lead  directly  to  chaotic  flows  [4]. 


Tig.  4.  -  Time  evolution  of  the  Fourier  component  6,,  of  the  temperature*  field  in 
two-dimensional  Biinard  convection  at  i?,  =  120Jf„.  I’r  =  20  [3],  The  numljcrs 
lalcelliug  each  plot  give  the  wave  number  cut-off  used  to  derive  the  Galerkin  approxima¬ 
tion  to  the  Botissinesn  equations.  Thus  1  ;'.2  gives  the  Lorentz  equations,  while  tlie 
higher-order  models  are  higher-order  Galerkin  approximations.  Observe  that  as  the 
re.solution  inereases  the  chaos  disappears. 
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2'J.5.  S|iniioiis  ( II  ii  iiicri  ra  1 )  t  ii  ib  n  1<‘ ii  »••• .  Orititv  it  a/.  >linw 
Ibat.  M  liili'  luw -order  dyiiatiiii‘:il  >yAlt“nis  dcriM-tl  Ity  (JaJfrJiiii  aj)|)ro.\ijiialioJi 
to  the  t  \vo-djiiieii>ion;!i  lioii.ssine.si)  eiiuatioli.s  may  exliiliit  eliaotie  xiliit ioii.s, 
this  ehaos  ty|iieally  disa|)i>ears  as  llie  (limeiLsioii  of  tli<“  projection  sjtr.ee  increases 
fsc<'  lie.  11.  Similarly,  it  was  sliotMi  by  OitszAu  iiml  Kklls  (5J  that  iimlcr- 
resolvcd  numerical  ealculatioiis  of  transitional  ])lanar  shear  flows  may  be  spu¬ 
riously  chaotic.  Umler-rcsolved  computations  do  not  have  de;rrees  of  freedom 
associated  with  small  spc.tial  scales  available  to  act  as  an  eddy  viscosity  on 
well-res(dved  larj'e  scales. 

2'3.  Computer  simulations  of  turbulence.  -  In  this  subsection,  I  shall  give 
three  examples  of  numerical  simulations  of  turbulent  flows.  The  first  two 
examples,  turbulent  channel  flow  and  the  simulation  of  a  turbulent  spot,  are 
of  the  natme  of  numerical  experiments  in  which  the  numericist  uses  the  computer 
in  much  the  same  way  as  the  experimentalist  uses  the  laboratory,  namely  as  a 
source  of  data  about  flows  in  a  controlled  environment.  The  final  example, 
the  Taylor-Green  vortex,  is  an  example  in  which  the  computer  is  being  used 
to  try  to  uncover  fundamental  physical  laws  of  turbulence. 

2‘3.1.  Turbulent  channel  flow.  Tmbulent  channel  flows  have  been 
simulated  numerically  three  ways:  a)  large-eddy  simulation  with  a  subgrid 
scale  tm  bulence  closure  for  eddies  outside  the  wall  layer  and  a  heuristic  boundary 
condition  applied  at  the  edge  of  the  viscous  sublayer  by  Dk.4KDorff  [11]  and 
Schumann  [12],  6)  large-eddy  simulation  with  a  subgrid  scale  turbulence  clo¬ 
sure  applied  to  eddies  of  all  scales  including  those  in  the  wall  layer  by  MoiN 
and  Kqi  [13]  and  c)  full  numerical  solution  of  the  Kavier-Stokeg  equations  by 
Orszag  and  Patera  [14].  The  really  crucial  differences  are,  as  we  again  note 
in  subsect.  2‘4  below,  between  a)  and  6)-c).  Simulations  of  type  a)  have  much 
smaller  computational  requirements  at  a  given  Reynolds  number  R  than  either 
of  types  h)  or  c),  the  latter  requiring  asymptotically  similar  computational  work 
at  large  R.  The  deficiency  of  simulations  of  type  a)  is  that  they  require  model¬ 
ling  of  wall  layer  effects  in  terms  of  an  over-simplified  boundary  condition;  the 
deficiencies  of  types  b)  and  c)  are  that,  with  currently  available  computer  reso¬ 
lution  (say  64x64x65  on  a  Cray-1  computer),  Reynolds  numbers  are  limited 
to  about  10  000  (type  6))  or  5000  (type  c)).  For  simulations  of  types  b)  or  c), 
the  computational  work  scales  as  R*,  so  future  increases  in  computer  power 
do  little  to  increase  the  effective  Reynolds  number  of  the  computations. 

Xevei  theles.s,  it  is  possible  to  achieve  interesting  results  with  full  numerical 
solutions  of  the  Xavier-Stokes  equations.  In  fig.  5,  we  plot  the  mean  velocity 
profile  found  in  the  channel  flow  computations  of  Orszag  and  Patera  [14]. 
The  fit  to  a  logarithmic  wall  layer  velocity  profile  is  only  marginal,  but  the 
resulting  von  Karman  constant  0.45  is  witliin  experimental  bounds,  so  this 
calculation  does  give  the  first  computation  of  a  wall  layer  from  the  basic  prin- 
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Fig.  5.  -  Mean  turbulent  profile  obtained  by  full  numerical  siinulation  of  plane  Poi- 
aeuille  flow  at  JJ  =  5000  using  a  64  x  64  x  65  spectral  simulation.  Xote  the  viscoiu 
sublayer,  buffer  region  and  logarithmic  layer  of  8-^9  data  points  (from  [14]). 
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Fig.  0.  -  A  plot  o!  the  turbulent  wall  pressure  spectrum  as  a  function  of  frequency 
(from  [15]). 


ciples  of  fluid  dynamics.  Another  more  reecnt  result  from  computations  of 
this  type  is  given  in  fig,  6,  in  which  we  plot  the  wall  pressure  spectrum  in  a 
moderate-resolution  (32  X.32  x33)  run  compared  with  available  experimental 
data  (see  [15]).  Despite  the  moderately  low  Kcynolds  number  {R  =  5000)  of 
tlie  simulation,  agreement  i.s  achieved  because  flow  features  that  do  not  depend 
explicitly  on  the  boundary  wall  layer  structure  tend  to  be  Reynolds  number 
indciieiulcnt. 
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2  •>.-.  'r  11  r  ti  11 1  «■  n  t  >))<»(.  'rinTi’  lias  hci'ii  iiiin  h  nii'iit  inIcirNi  In  lie 
cMililt  inn  of  liualixfil  «s|iols»ili  I  iirliulrllt  (lows  (.sro  |  Mi  |).  TlirliiNt  liiillii'Mra  I 
.>iliiiil:'lioli  of  a  ( iirtiiili'iil  .spot  was  rf)to)  ioil  liy  JjaiVAIllt  i  1  7  |.  wlio  ii.vcd  thici- 
iliniclisioiial  vortex  tilalin-iil  teclmiiiiirs  to  l•(llllJllll»■  the  (invisiiil)  How.  .A[ric 
recently,  we  have  be;;un  a  stuily  of  .spots  iisiiie  full  iiimieriial  selutioiis  of  the 
Xavier-Stokes  eipiatioiis  at  moderate  Keynohls  miiiibers  [IS].  The  latter  simii- 
lafioiis  are  iierfcrtneil  by  forein<;  the  initial  How  nsin^:  a  localized  fun  e  to  di  i\e 
a  jet  of  fluid  vertically,  then  allowiii};  the  disliirbai:ce  to  evolve  naturally. 
In  H<;.  7  and  S,  we  plot  contours  of  muxiimim  xerlica.l-c  velocity  in  the  (.',  //) 
and  yr,  plane.s  at  various  times  of  evolution  of  plane  roiseuille  How.  The 
charact«‘r  of  this  siwt  evolution  is  similar  to  that  observed  experimentally; 
the  sjKit  seems  to  spread  in  the  spanwise  direction  by  *  fransver.-e  contami¬ 
nation  *,  in  .a/rreement  with  the  dye  injection  experiments  of  Gad-el-Hak 
et  al.  [19];  the  ftreatest  turbulent  activity  is  near  the  edfres  of  the  .sjiot,  the 
« .spreading »  angle  of  the  spot  relative  to  it.s  source  is  about  1(»^,  in  agreement 
with  the  channel  flow  experiments  of  Carlson  et  al.  [20];  the  vertical  structure 
of  the  spot  i.s  in  qualitative  agreement  with  that  observed  experiment ..lly. 
Further  numerical  experiments  are  under  way  that  should  elucidate  deti.ils  of 
the  flow  ill  spots  and  tlie  surrounding  fluid. 
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Fig.  7.  -  Contours  of  the  tnaximum  s-velocity  in  the  (j:,  y)-plane  at  t  =  12,  18  after 
initializing  a  turbulent-spot  computation  by  imposed  vertical  forcing.  These  computa¬ 
tions  are  performed  using  a  spectral  code  with  128  (jr)  x32(y)  x32(z)  resolution.  Fourier 
series  are  used  in  x  and  y;  Chebyshev  polynomial  expansions  are  used  in  r.  Here 
li  =  GOOO. 
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Fig.  8.  -  Same  as  flg.  7,  except  x-z  contours  of  inax,lc,j. 

2‘3.3.  Taylor-Green  vortex.  In  order  to  gain  understanding  of  the 
basic  jdivsics  of  the  generation  of  small-scah*  turbulent  flow  features,  a  nice 
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III!  licl  IH'olilcIll  In  till-  T.l.vlo|'-<il  <■<•11  ('r<i)  MlltcV  llllW  | 1  .  'J  Ll  J.  Ill'll'  till'  t!llU  i' 
tliiit  uliiili  (l<-\<‘l()jis  ill  tiiiii-  friiiii  initial  ((iinlil  i(iii>  tlial  ((iii>i.'t  nl  i-Ncilal  inn 
in  liii'icalls  a  sin.u'li'  KnuriiT  ninili'.  Ui'iaiisi'  nf  ininlin<‘ar  interact imi,  the  llnw 
liei  niiies  .strnii;rly  three  iliiiK  iisiniial  and  develnjis  <‘X(  italinn  at  all  sjiatial 
seali-'.  The  T(1  vnvtex  I'.as  lieen  nseil  to  study  such  fuiulaineutal  iiuestinns  as 
the  i".ihaMeeini'nt  <if  vortieity  liy  Mirtex  line  stretehinn,  the  a]iiir<iaeh  to  isoti’oiiy 
of  the  'inall  scales.  ])().ssil)l<*  siniriilar  Ixdiavior  of  the  Euler  equations,  forniation 
of  an  inertial  ranj^e  and  analysis  of  the  ;:eoni“try  :ind  interniitteney  of  hi^di- 
A’ortieity  recioiis.  The  TG  tloiv  is  advantaireous  for  these  studies  because  its 
siM*cial  symtnetry  has  allowed  the  development  of  numerical  aljjorithms  that 
are  a  factor  64  more  efficient  in  both  memory  and  storaffc  than  conventional 
periodic-geometry  spectral  methods.  For  a  three-dimen.sional  flow,  this  faitor 
64  translates  into  a  factor  4  increased  range  of  spatial  scale- — it  is  now  pcs.'iblc 
to  compute  the  TG  vortex  flow  with  512x512x512  Fourier  modes  for  each 
velocity  component  on  the  Cray-1  computer  (or  more  than  4-16'  effective 
degi'ces  of  freedom  f). 


rig.  U.  ~  A  plot  of  the  dMtrihutioii  of  large-vorl icily  regions  m  the  TG  vortex  Ihnv 
ns  a  fiiiirtioii  of  time  I  and  ilistaiire  </  away  from  the  side-walls  of  the  impciineahle 
ciihe  in  Mhich  the  How  takes  )>laoe.  tllisevve  how  vortieity  explodes  in  towards  the 
<’eiiti-r  of  the  ciihe  l>etH<-<'ii  /  4  .iiiil  I  S  /friiiii 
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(tile  nl'  dll'  Iiiiiic  CM  Miti"  I<  'ill1s  In  ^•||n•l:;|•  hi  in  i m  slinlii  >  <,|  d.c  T'i  lliiw 
i.'  I  lie  .■<l!”';;<->di»li  I  lint  \  l^<•(,•^it^  lii;i  y  jil.iV  .ill  <'.->ilil  I;' I  i  rlc  Ill  die  iii'\  i  .(ili:ii  I  nt 
.smilll-sciilc  tUI'lillh-lK'f,  liol  jli>t  liitili;;  ns  il  sink  nf  dirlnilclit  klllctlr 
Indootl,  \vn  lind  lliat  tlic  lU-M-lnpnu'iit  nf  die  I urlnilcnt  llnw  scviiis  to  iciiuiin 
viscnsif  y  In  iiidin c  iiistaltilidc.s  of  Mirdci  l  sImicIiim  s  in  liii  )i  tlic  inil  ial  )iir;;c- 
scaU“  uonturliulcnt  vordcity  uiidt  i'onfs  an  cxjilrsivn  rodisf  i  dnitinii  in  spate 
(set-  ti;:.  9).  Those  vi.sfcsity-indiifoil  instabililios  aie  jirnhahly  oliot  five  bocaiiM* 
visco.sity  allows  vortex  line  recoiineff ions  jirtdiibiteil  in  inviscid  flow.  Sinhlar 
diiTusioiial  instabilities  liave  now  iK-en  shown  to  be  res]K)nsible  for  the  ■'eiieration 
of  small-seale  .struetures  in  two-dimensional  maenetohydrt.dynamie  [23]  and 
kinetic  [24]  turbulence.  Further  study  of  viscosity-iiulin  ed  instabilities  should 
clarify  tin*  tlmelopnient  of  intermittent  flow  strui  tnies  in  turbuleme. 

2  4.  Suhgrid  »calt  (urbvience  clvuurei).  -  I’erhaps  the  most  distin;;uishiiij;; 
characttristic  of  hifih-Reynolds-number  turbulent  Hows  is  their  h.rjre  raiipe  of 
excited  space  and  time  scales.  In  lK)mo;;eneous  turbulence,  dissipation  scale 
eddies  are  of  order  It'  times  smaller  than  encrfry-containinir  eddies.  Includin;; 
the  effect  of  this  ruiiK^  of  spatial  scales  on  the  allowable  time  step  in  a  numerical 
solution  of  the  Xavier-Stokes  equations  gives  the  estitnate  that  order  E*  oper¬ 
ations  are  required  to  simulate  a  turbulent  flow.  This  is  the  reason  for  interest 
in  the  large-eddy  simulation  method  in  which  excitations  on  scales  smaller 
than  those  resolvable  numerically  are  modelled,  usually  by  an  eddy  viscosity 
coefficient  (see  [11, 12]).  The  basic  action  of  an  eddy  viscosity  on  large  eddies 
is  reasonable,  although  it  cannot  reproduce  the  random  character  of  the  action 
of  small-scale  eddies.  However,  in  order  to  model  projierly  wall  turbulence, 
it  is  necessary  to  extend  the  subgrid  modelling  ideas  of  Deardorff  and  Schumann 
and  treat  the  turbulence  all  the  way  up  to  a  rigid  wall,  as  in  recent  work  by 
Moin  and  Koi  [I'l],  Unf«  rtunately,  in  order  for  Moix  and  Km  to  resolve 
motions  down  to  the  scale  of  turbulent  bursts,  which  is  necessary  in  order  to 
capture  the  meehaui.>m  producing  the  turbulence,  the  work  restriction  0(li*) 
remains.  Thus  the  Reynolds-number  restrictions  are  similar  for  large-eddy  and 
full  numerical  solutions  of  the  Xavier->Stoke.s  equations  that  attempt  to  inte¬ 
grate  all  the  way  through  the  wall  layer  region. 

In  ref  at  work,  Yakhot  and  Orszag  [25]  have  used  dynamic  renormali¬ 
zation  group  (RXG)  methods  to  treat  wall-bounded  turbulence.  The  idea  of  the 
infra-red  RXG  method  is  to  use  perturbation  metheds  based  on  the  direct - 
interaction  approximation  [26]  to  eliminate  all  small  spatial  scales  up  to  the 
resolvable  grid  scale  from  the  Xavier-Stokes  equation.s.  Tins  is  done  iierfiirh;!- 
tively  by  eliminating  narrow  baml.s  of  wave  vectors  from  the  dynamics  (kc 
fig.  10),  renormalizing  the  resulting  rediued  dynamical  equation  to  l.ave  the 
form  of  the  Xavier-Stokes  equation  with  tnoditied  viscosity  and  ramlcin  forcing 
ternis,  and  then  repeating  the  process  iteralivelr-  until  all  the  reijuircd  small 
scales  are  removed.  The  resulting  dyinnnhal  etniations  involvf'  a  moditied 
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oililv  \i>(i.>ily  amt  :i  ramlinii  fcnc.  Imlli  iiiiliicoil  liy  rein  riiiali/.aliiui.  'I'lie 
eilily  viNi  ii-iily  is  iiiinlitieil  friiiii  (lie  Siiia;;(iiiiisky  vixesity  umiI  liy  ] t);Ai;iiiii:i  r. 
SCUOIANN  anil  ^loiN  amt  Kim  in  (lie  wall  regions  in  wliieli  there  is  iiiterferein  e 
lielween  the  eihly  anil  imileeiilar  A  iscnsilies.  Tliis  interf<  ren(  <■  elTeel  is  the  ke\' 
to  olitainiii;,'  ;i  ftiitlifiil  re|ire.sent:;tion  nf  the  wall  region.  Also.  Ilie  iiniineil 


0  /r,  k  /le»p[-i.]  A 

I'i;:.  10.  -  A  M'ticniatic  representation  of  the  uiodiil  sfrneture  of  (be  dyiiaiuic  renormal¬ 
ization  group.  Here  represents  wave  numbers  witbiii  tbe  energy-coiitaiiiing  range, 
while  .1  gives  the  bigh-wave-number  (viscous)  cut-off.  Modes  in  tbe  batebed  band 
are  removed  at  each  step  of  tbe  RXG  procedure. 


random  forte  is  large  in  the  buffer  layer  between  tbe  viseou-s  sublayer  and  the 
logaritlimic  layer,  giving  a  turbulence  source  in  this  region.  Further  work  is 
now  under  way  applying  these  KXG-based  closures  to  both  large-eddy  simula¬ 
tions  of  turbulent  shear  flows  and  to  the  deriv.ition  of  new  classes  of  turbulence 
transport  (Iteynolds  averaged)  equations  that  should  be  useful  in  engineering 
applications. 


3.  -  Conclusion. 

I  have  reviewed  several  areas  of  activity  in  the  numeiieal  simulation  of 
transition  and  turbulence  in  which  I  have  been  intimately  invohed  recently. 
In  this  short  space,  it  has  not  been  possible  to  do  justice  to  all  of  the  large  number 
of  researchers  involved  in  these  fields;  the  references  do  a  more  complete  job  of 
surveying  the  literature.  The  principal  conclusions  from  our  studies  are: 

i)  Xumerical  methods  now  provide  essential  information  complemontarA* 
to  that  available  from  experiment  and  mathematical  analysis. 

ii)  Comput.ational  fluid  meclianic.s  ha.s  now  matured,  so  that  tlure  are 
techniques  that  can  be  reliably  applied  to  the  most  difficult  of  fluid-mcchauical 
problems.  In  contrast  to  10  years  ago,  it  is  no  longer  mainly  a  question  of 
how  to  compute  a  complicated  flow,  rather,  now,  it  is  a  question  of  which  flow 
to  compute  in  order  to  extract  tlie  most  useful  information. 

iii)  It  i.s  cruci.il,  e.speeially  in  our  studies  of  transition:;!  flows,  that  we 
have  used  spectral  numerical  methods  (sec  sect.  1  above).  Spectral  methods 
are  so  accurate  for  these  problems  that  we  can  confidenth'  conclude  that  jiroji- 
erly  tested  numerical  result.s  arc  true  fluid-mecluinical  results.  In  contrast  to 
finite-dilTereiu  e  or  fmite-elenient  methods  in  which  an  increase  in  spatial  resolu¬ 
tion  by  a  factor  2  leads  to  an  error  dei  isuii-e  by  a  factor  I  or  IS  or  so,  w  itli  sjx'ctral 


nu'tliiiils  :i  fiictiii'  -  ili<  r<-:'sc  in  r«  >oliit imi  f  \  |pir:tlly  <l<(  s  (luMTicr  liy 
(ii'di  rs  (if  utl<‘.  Tliis  •.icr  iinitt'  \ t  l  ilicnl inn  of  icMiits.  For  cximiidc, 

in  ifcoiit  .sliiili<s  of  tiMii.'ilioii  in  Conctlo  flow,  M.\itrrs  <t  id.  j:,'!]  and 

M Ai:ci's  [L'S]  luivo  luM-n  aide  to  indiiovo  at  least  tlim‘-de(  imal-j)lace  ajna'eiiieiit 
witli  e\i)oriinent  on  wave  sikmhIs.  The  confidence  in  these  results  ha-  pciinittc<l 
new  analytical  insights  into  the  charai  ter  of  the  ons*‘t  of  wavy  instahilit i<-s  of 
Ta>  l<M'  A'ortices  in  Couette  fi<)W  I'Jtt]. 

iv)  Xew  venerations  of  hi;:ver  and  faster  coinputers  can  most  prolitahly 
be  Used  to  extend  the  raiifre  of  applicatioji  of  computational  Hiiid  dynamics. 
Transition  and  turbuleiue  probh-ms  in  complex  geometries  with  complex 
])hysics,  like  multiphase  flows,  will  surely  be  the  subject  of  studies  in  the  near 
future. 
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